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ABSTRACT 

Most of the well known algorithms for watermarking of digital images involve transformation of the image data 
to Fourier or singular vector space. In this paper, we introduce watermarking in Hilbert transform domain for 
digital media. Generally, if the image is a matrix of order m by n, then the transformed space is also an image 
of the same order. However, with Hilbert transforms, the transformed space is of order 2m by 2n. This allows 
for more latitude in storing the watermark in the host image. Based on this idea, we propose an algorithm 
for embedding and extracting watermark in a host image and analytically obtain a parameter related to this 
procedure. Using extensive simulations, we show that the algorithm performs well even if the host image is 
corrupted by various attacks. 
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I. INTRODUCTION 

A large amount of accessible information today is available in one or the other multimedia formats. While this serves the 
purpose of easier dissemination of data, this also makes it vulnerable for misappropriation and misuse. Hence, it is all the more 
important to protect intellectual property rights of the content available in various digital formats. Digital watermarking IS a 
popular method by which the owners of the data, in any multimedia format, can embed their logo, trademark or some proprietary 
information |j3l |4t] in a way that can either be visible or invisible to a general user. This information can be later retrieved for 
verification purposes or in case of conflicting claims on the ownership of data iH-Q]- 

Thus, when applied to the case of digital images, digital watermarking technique consists of (i) an algorithm to embed a 
watermark image on a host image and (ii) an algorithm to retrieve the embedded watermark with least distortion. Ideally, we 
would expect that the algorithms be robust against any manipulation of the original data and it should also be designed to render 
any illegal retrieval of watermark a futile exercise. The field of digital watermarking has been the focus of research attention for 
more than a decade now. This is partly due to the proliferation of multimedia formats as well as new tools, both commercial 
and open source, to manipulate them. In this paper, we focus on the invisible watermarking of digital images using the Hilbert 
transform technique. 

Watermarking techniques for digital images can be broadly classified into two categories, namely, the spatial domain tech- 
niques and transform domain techniques depending on which domain the watermark is embedded. Typically, in spatial domain 
techniques the watermark is embedded in those part of the data that do not distort the host image in any significant way. For 
instance, some of the well-known spatial domain techniques are least significant substitution |8, 9] and the correlation based 
approach i lolfTTll . In least significant substitution technique, the watermark is embedded by replacing the least significant bits 
of the image data with the bits of the watermark data. There are many variants of this technique. In correlation based approach 
the watermark is converted to a pseudo-random noise (PN) sequence which is then weighted and added to the host image with a 
gain factor. For detection, the watermarked image is correlated with the watermark image. In the transform domain techniques, 
the watermark is embedded in those parts of the transformed host image which do not distort the image significantly. 

One of the earliest transform domain techniques is the one based on discrete cosine transform (DCT) [12-.- 16.1. In DCT, 
the image is decomposed in terms of various frequency bands and watermarks are embedded in the middle frequency bands 
which are not significant for the host image. Further, image transformations do not affect the watermark placed in those bands. 
DCT based methods are generally robust, particularly against JPEG and MPEG compression. The techniques based on wavelet 
decomposition are similar in spirit to DCT with the additional feature that the multi-resolution character of the wavelets allows 
graded information to be stored at various resolutions. For instance, in 1 17] wavelet coefficients of the image and the watermark 
at different levels of resolution are added together within the constraint of the so-called human-visual model. A review of 
wavelets based techniques is available in |[l8-22]. There is yet another method of digital watermarking based on singular value 
decomposition (SVD) techniques ll23l - [32ll . In contrast to DCT and wavelets based techniques, the advantage of singular value 
decomposition based methods is that they provide a transform space that is tailor made for the given image data matrix. Both in 
the DCT and wavelets, the basis for the transform space is a fixed set of functions. In the SVD, it must be calculated from the 
given data and the singular vectors so calculated form an optimal basis for the image matrix in the least square sense. It is worth 
mentioning that some authors have resorted to hybrid techniques i.e., algorithms based simultaneously on different domains to 
improve the watermarking results lf33l - l36ll . 

In this present work, we propose a new scheme for watermarking digital images using the Hilbert transform. The analytic 
signal s{t) associated to a signal s{t) is, 

i{t) = s{t) + isH{t) (1) 

where snit) is the Hilbert transform of s(t). Clearly, .f(f) can be written in phase-amplitude form as A{t)e'^^'\ If the signal 
changes sufficiently slowly, then the phase of the analytic signal is negligible. Typically, most images have slowly varying pixel 
values except at the edges. In such a scenario, we can expect the phase to be negligible most of the time and the matrix of phase 
values will be sparse. Hence a good amount of information can be embedded in the phase of the analytic signal associated to the 
image. Since only the phase of the analytic signal is proposed to be used for embedding the watermark, it is likely to be highly 
imperceptible for visual perception. This is one of the key requirements of ideal watermarking algorithms. In addition, this also 
provides a large space for embedding the watermark which is useful for creating redundant watermark distributed throughout the 
Hilbert transformed space. This makes the algorithm more robust against attacks. This is the main idea underlying the present 
scheme and to the best of our knowledge the Hilbert transform has not been used for watermarking purposes before. 

Further, our proposal is a non-blind scheme, which implies that to recover the watermark from the watermarked image, we 
require both the original as well as the watermark image. This is not necessarily a restrictive requirement as there are many 
situations in which this scenario is valid, such as in ownership litigations. In many such cases, the owner keeps a copy of both 
the original host image and the watermark. In addition, non-blind scheme also makes detection of watermark in an arbitrary 
image difficult if the watermark image is not available. 

The rest of the paper is organized as follows. In the next section, we introduce the proposed method of digital watermarking. 
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In Sec. Ill, we provide the formula to optimize the scaling factor Section IV describes the results of numerical simulation, while 
Sec. V describes robustness of the algorithm. Section VI gives a summary and concluding remarks. 



n. WATERMARKING USING THE HILBERT TRANSFORM 

Given any arbitrary signal s{t), as a function of time f , we can construct an analytic signal of the form 

s{t)^sit) + isHit)^A{t)e''^^'^ (2) 
where is the Hilbert transform of s{t). It is defined in terms of the Cauchy principal value (P.V) of the integral 

SH{t) - -P.V r ^ dr (3) 

provided this integral exists. In general, Hilbert transform has a wide range of applications, in particular, in the area of signal 
processing iIjtIi . In signal processing, the Hilbert transform zh of a discrete signal z is defined as the output of a linear filter with 
the frequency response H{(£)) given by 

{/, -n <(o<0 
Q<(0<K (4) 
0, (o = -n,Q,n 

In the context of this work, let z denote a signal vector of size n x 1 as a function of position (in arbitrary units). For instance, z 
would represent one column of pixel values of an image matrix Z of order nxm. Then, the analytic signal z associated with z is 
z + izH, a complex valued signal with phase and amplitude. We intend to embed the watermark image in the phase component. 
We denote the amplitude and phase of this analytic signal by n x 1 vectors and 9^ respectively. The inverse Hilbert transform 
consists of taking the real part of this analytic signal. Hence we have, 

z = rz*cos0z, (5) 

where the * stands for element-wise multiplication and is not to be confused with usual matrix multiplication. Further, for any 
vector 9z, the operation cos 9z represents the cosine of every element of the vector. 

Notice that the Hilbert transforms in two dimensions is not uniquely defined. Hence, in the watermarking algorithm presented 
below, we will treat the image matrix Z of size nxm as n vectors and apply Hilbert transforms to each of the n vectors. From 
a computational point of view, for the Hilbert transform of any image Z, we are required to compute n independent Hilbert 
transforms each with m elements. However, for ease of notation, we present the algorithm using scalar elements instead of 
in vector notation. Now we will apply this to define Hilbert transform of an image matrix Z of size n x m with elements 
Zjj,i — 1,2,3, .. .n,j — 1,2,3, .. .m. Hence Eq.® will still continue to hold in the form 

Zij — aijcos&ij, i — 1,2,3, .. .n, j = 1,2,3, .. .m. (6) 

In this, Uij are the amplitudes and 0,y are the phases obtained from vector-wise Hilbert transform applied on Z. 

A. Algorithm for embedding the watermarli 

In this section we describe the algorithm to embed a watermark image into another gray scale image of the same size using the 
Hilbert transform. Let the host image be represented by Z with elements Z,y and W denote the matrix of the watermark image 
with elements Wij. 

First, we perform the Hilbert transform of the original image Z and obtain the relation 

Z,-^' =fl,;COS0i^'. (7) 

Next, we do the same for the watermark image W to get 

Wij = b,jCos<Pij, (8) 

where bij and <I>,y are the phase and amplitude. 
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Now, we add the scaled amplitude bij of the watermark to the phase &ij of the original image to get 

eij = &ij+Xbij, (9) 

where A is the scaling factor For typical images, the order of magnitude of bi/s is much larger than Qifs for all i,j and hence 
A will have to be much smaller than unity in order to compensate for this difference in order of magnitudes. 
Finally, we get the watermarked image as 

ZYj = ai j cos Oij. (10) 

Thus equations (l7b-(fT0li constitute the algorithm for watermarking using the Hilbert transform applied column-wise to an image 
matrix. We remark that this algorithm is motivated by the fact that most of the information about the image is encapsulated 
within the amplitude. The phases 0,y, for all / and j, contain very little information and can be thought of as a sparse matrix. 
Hence, we can store most of the information about the watermark without causing too much distortion in the original image by 
following the above strategy. 



B. Algorithm for extracting the watermark 

Given the watermarked image Z™, we can extract a (possibly corrupted) watermark if we have access to 0,y,<I>,j, fl,j, for all 
/ and j, and the value of X. This information is most easily available if one has access to the original image as well as the 
watermark image. As pointed out in a previous work of the first and the third author (s^l, this is not a particularly restrictive 
assumption. 

The extraction algorithm is just the reversal of the embedding algorithm given in the previous subsection. Starting from Eq. 
([Tol l we divide both sides by aij and use the inverse cosine function to recover Oij. By substituting for dij using Eq. (|9]l we get 
for the amplitude of the watermarked image 

^ cos-' (Z^. / at 

bij^ ■ (11) 

Finally using Eq. (O, the extracted watermark image can be constructed as follows, 

Wij^bijCos{<Pij). (12) 
Thus, Eq. ( fT2] i. along with Eq. ( fTTT i constitute the watermark extraction algorithm. 



III. OPTIMIZATION OF THE SCALING FACTOR 

The choice of the value of the scaling factor A plays an important role in our watermarking algorithm. If A is chosen too 
small, then the quality of embedding is good but that of the extracted watermark is poor. On the other hand, if A is too large, 
then extraction works well but embedding suffers. Hence, A needs to be chosen optimally to achieve a balance between these 
extremes. In general, this problem is a non-linear optimization problem and we shall describe an iterative method to produce a 
solution to it. 

We shall use the mean square error (MSE) as a measure of the quality of embedding or extraction, and attempt to minimize 
the sum of the MSE from the embedding and extraction steps. First define the function /(A) as 



/(A) = ^^ (Z--Z„-)2 + (W,;-W„)2 



(13) 



With this notation, we need to find the value of A that minimizes /(A). Note that both Z"" and W depend on A. We also note 
that, in practice, converting to an image introduces a truncation error £, which, though negligible in the embedding step, 
affects the extraction step significantly. 

Now, we use Eqs. (|9]l and (flOt adjusted for the truncation error e to get the following 

m n m n 

^ i:(Z;; -Z„)2 = ^ ^ [a^jCos{@,j + lb,j) 

,7=1 i=l 7=1 1=1 

+ e-aijCos@ij\^ (14) 

m n 
7=1 (=1 
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where the terms involving e and higher powers of X are ignored. 
Similarly we use Eqs. (fTTT i and (fTSl i adjusted for e to get 

m n 



/n n 

;=i'=i 



cos ' [ COS {<dij + Xbij) + e /aij\ — @i 



X cos — bij cos <Pij 



EE 

j=U=l 



a(A)cos4>,-^ 



where a (A ) is given by 



a(A) = cos '(cos(0,-j + A/7,;) + e/ai;)-(©,7 + A/7,;)- 



(15) 



(16) 



With this preparation we return to the minimization of /(A). Our approach is to start with any initial value A = Aq and produce 
successive values Ai, A2, . . ., each depending on the previous value, which converge to the desired minimum of /(A). 

The method to produce A^+i from X(; at the £-th step is as follows. First we use Eqs. (fT4l i - ( fTSl l to approximate /(A) by // (A) 

as, 



2i -2 



fti^) = E E {aiibijSm@ij)'^X^ + {a (Xe) cos (i>ij)^X 

j=ii=i 



(17) 



Note that /f (A) depend on X(. Next, we choose X(^i to be the value at which /f (A) attains its minimum. After a straightforward 
computation, we have the following, 

V\ 1/4 



^e+i — 



Y^{a{Xe)cos^ij) 



(18) 



L(fl,,A7sin0,y)2 

Thus Eq. ( fTSl l can be iterated for £ = 0,1,... starting from any initial value Aq till the desired level of convergence is achieved. 
In practice, this algorithm seems to have fast convergence. For example, for our test images it converges within £ = 2 steps to 
A = 0.0018 (upto 4 decimal places) for a wide range of initial values. 



IV. NUMERICAL SIMULATION 



In this section we will apply our embedding and extraction algorithm to the host image shown in Fig. |l(a)| and the water- 
mark image in Fig. |l(b)| Both these images have dimensions 512 x 512. As pointed out earlier, the quality of the water- 
marked image improves and that of the extracted image deteriorates as the scaling parameter A — > 0. This is borne out by 
the simulation results and can be clearly seen from the watermarked images shown in Figs. |2(a)[ |2(c)l |2(e)| 2(g)[ a nd|2(i)|for 



A = 0.1,0.01,0.0018,0.001,0.0001 respectively. The corresponding extracted watermark images are shown in Figs. 2(b)[|2(d)"] 



|2(f)[|2(h)1 and 2(j) The simulations clearly demonstrate that the Hilbert transform based algorithm proposed in Eq. (I7b-(ll0b aTid 



Eq. (fTTll.(fT2l) produces results whose visual quality is good and acceptable as a practical tool for watermarking digital images. 

These results are quantified by means of peak signal to noise ratio (PSNR) and the root-mean-square error (RMSE) of the 
corresponding images. PSNR is used to evaluate the perceptual distortion of the proposed scheme. PSNR and RMSE is computed 
for the image difference matrix Ajj = lij — Uj for several values of A . In our context, there are two possibilities, (a) Ijj are the 
elements of host image and /,,■ are the elements of watermarked image, and (b) /,y are the elements of watermark image and /,y 
are the elements of extracted watermark. Then RMSE is defined as 



/I m n 



We also define the PSNR, measured in decibels, as 

_2 1^2 



p=101ogio(maxz;;^./e^), (20) 

where max zt.j represents the maximum value of a matrix whose elements are zt.j- Figure|3]depicts the combined RMSE versus 
lambda graph of the extracted and watermarked images. Notice that A = 0.0018 corresponds to a minima as predicted by the 
analysis in the previous section. 
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FIG. 1. |(a)| Host image and |(b)| watermark image of size 512 x 512. 



V. ROBUSTNESS OF THE ALGORITHM 

An important property of the watermarking algorithms is that they should be robust against various possible attacks. In this 
section, we will first embed the watermark in the host. The watermarked image will then be subjected to some attacks. After the 
attack has been carried out, the extraction algorithm is applied. We will subject our algorithm to the following major attacks (i) 
robustness against Additive Noise (ii) Cropping (iii) Gaussian Noise (iv) JPEG and JPEG 2000 compression (v) Median Filter 
(vi) Rotation (vii) Gamma Correction (viii) Intensity Adjustment (ix) Gaussian Blur (x) Contrast Enhancement (xii) Dilation 
(xiii) Scaling. The full list of attacks and their result is given in Table. (HI. Before we proceed, at the very outset let us mention 
that all the images are taken to be of 512 x 512 and the scaling factor is chosen to be 0.0018 (optimal value). In particular, all 
the attacks are done on Fig. |2(e)| 



A. Additive Noise 

Our first attack will be the addition of uniformly distributed noise to Fig. |2(e)| resulting in the corrupted image shown in Fig. 
|4(a)| Figure [4(b)| shows the result of extracting the watermark. As can be noticed visually and also from the p and e values, the 
quality of the extracted watermark is quite good. 



B. Cropping 

Cropping is the process of removing a portion of the image. In the present case we will remove a portion of the watermarked 
image and check if our extraction algorithm is able to extract the watermark. Figures [5(a)| and [5(b)| show the cropped watermarked 
images and the extracted watermark respectively. Though there is some residual effect of the host in the extracted image, it is 
worth noting that it appears only in the cropped portion of the image. The extracted image is unaffected by this residual effect of 
the host. We would also like to mention that we have carried out cropping for various percentages (10%, 20%, 30%, 40%, 50%) 
of Fig. |2(e)| and the extracted image is good. 



C. Gaussian Noise 

When Gaussian distributed noise is added to the watermarked image we obtain Fig. |6(a)| The extracted image from this is 
shown in Fig. |6(b)| Once again our extraction algorithm performs well and yields the watermark. 



D. JPEG and JPEG 2000 Compression 

JPEG compression jssll is one of the popular standards for encoding of digital images. This is also an important test that 
any watermarking algorithm has to defend itself against. We subject Fig. |2(e)|to various percentages of JPEG compression and 
consequently extract the watermark from it. Figures [7(a)[ |7(c)[ |7(e)[ and |7(g)| are JPEG compressed at 20%, 40%, 60%, and 
80% respectively. The extracted watermarks from the aforementioned compressed images are shown in Figs. |7(b)[ |7(d)l |7(f)[ 
and |7(h)| respectively. From the PSNR values indicated in the table, it is clear that the quality of the images extracted as well as 



7 




FIG. 2. This panel shows the watermarked image on the left panel and the extracted watermark on the right panel for |(a)[[(b)l A = 0.1, |(c)[[(d)l 
1 = 0.01,[(e)][(2]A = 0.0018,[(g)][(h)]A = 0.001 and[(I)l[U)]A = 0.0001. 




FIG. 6. |(a)| Gaussian noise watermarked image and | (b) [ Extracted watermark. 
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FIG. 7. This panel shows the watermarked image on the left panel and the extracted watermark on the right panel for |(a)[[(b)] jPEG compressed 
20% watermarked image and Extracted watermark |(c)[|(d)| JPEG compressed 40% watermarked image and Extracted watermark |(e)[|(f)| JPEG 
compressed 60% watermarked image and Extracted watermark ) (g) [ [(h)] JPEG compressed 80% watermarked image and Extracted watermark. 



the embedded is better when the compression is less. This seems normal since JPEG is a lossy compression method, the higher 
the compression, the more is the loss in information of the image and vice versa. 

Now we focus on JPEG 2000 standard and show the results of attacks using JPEG 2000. This is a more recent standard 
for image compression and is based on wavelets and is much more efficient in reducing the size of images. Figure s [8(a)| and [8(b)| 
show, the compressed watermarked image and the result of the extraction from it respectively at a JPEG 2000 ratio valued at two. 
The PSNR and RMSE plot for the watermarked image for various JPEG 2000 compression ratios is given in Fig. |9(a)| Similarly 
the plot for the extracted watermark is given in Fig. |9(b)| Just like the JPEG case, the results show a deterioration in the quality 
of the watermarked as well as of the extracted watermark images as the compression ratio is increased. However, the watermark 
can still be obtained even for very high compression ratios which could be attributed to the robustness of the algorithm. 
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E. Median Filter 



Median filtering is used to remove outliers without reducing the sharpness of the image. It is similar to an averaging filter 
in which the value of the output pixel is the mean of the pixel values in the neighborhood of the corresponding input pixel. 
However, in the present method of median filtering, as one might have guessed, the value of an output pixel is determined by 
the median of the neighborhood pixels. The extracted watermark, Fig. |10(b)[ from the median filtered Fig. |10(a)| shows that the 
proposed scheme is robust to such an attack. 




FIG. 10. |(a)| Median filtered watermarked image and |(b)| Extracted watermark. 
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(a) (b) 
FIG. 1 1 . |(a)| 1 Degree rotated watermarked image and | (b) [Extracted watermark. 




(c) (d) 

FIG. 12. |(a)[|(b)|Gamma 0.95 watermarked image and Extracted watermark) (c) [ |(d) [Gamma 1 .09 watermarked image and Extracted watermark. 



F. Rotation 



Figure fTl (a)| shows the result of rotating the watermarked image. Fig. |2(e)[ through one degree in the anti-clockwise direction. 
After this has been done we crop the four corners of the rotated image in order to keep the same size as the original image. Now 
when the extraction algorithm is applied, the resulting extracted watermark is shown in Fig. |1 l(b)| Tests have been done for 
rotations by 2 and 3 degrees too. In these cases we are able to retrieve more than 50% of the watermark after extraction. 



G. Gamma Correction 



In general many output devices, say computer displays, have a nonlinear response to that of the input signal. Gamma correction 
is applied to compensate for this effect. Typically, the response follows a power law; . In the case of images, a gamma 
correction to it has the effect of compressing and expanding the pixel value. When 7 < 1 it is called compression and for 7 > 1 
it is called expansion. Since, 7 is related to the value of the pixel, which in turn directly translates into visual quality of images, 
a gamma correction leads to brighter or darker images depending on the value of 7. 

Figures |12(a)| and |12(b) show the watermarked image and the extracted watermark respectively after gamma compression 
with 7 — 0.95. Figures 1 2(c)| and 1 1 2(d)| show the watermarked and extracted images respectively after gamma expansion with 
7=1 .09. In both the cases, the extracted watermark is clearly visible and PSNR value is within the acceptable limits. 
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H. Intensity Adjustment 

Intensity adjustment is a technique of mapping an image's intensity values to a new range. In the present example, we have 
mapped the watermarked image's intensity to the range 0.1-1.0. As shown in Fig. |13(a)| and |13(b)| the watermark recovered 
from the intensity adjusted image does not suffer much distortion. 



I. Gaussian Blur 

The Gaussian blur filter, as implied by its very name, blurs objects. This filter creates an output image after applying a 
Gaussian weighted average of the input pixels around the location of each corresponding output pixel. Figure |14(a)| is the 
Gaussian blurred watermarked image and Fig. |14(b)| is the watermark extracted from it. 

Blurring can also be achieved by a simple low-pass filter Its mainly applied to average out the rapid changes in the intensity 
of the image. This is achieved by calculating the average of a pixel and all of its eight immediate neighbors. This average is then 
used instead of the original value of the pixel. This process is repeated for every pixel of the image. Since, this is very similar 
to the Gaussian blur attack, already presented, we will not give the resulting images of this experiment. Suffice to say that the 
extraction of the watermark is possible. 



J. Contrast Enhancement 

Contrast enhancement is the process by which the contrast of an image is modified. This can significantly improve the visual 
quality of the image. In this process the darker colors are made more dark and the lighter regions more light. This is achieved 
in the following way. We choose a particular cut off for the darker as well as the lighter regions and all the pixel values that are 
smaller and larger than this cut off are set to their mnimimum and maximum values respectively. For example, in the case of 
a greyscale image all pixel values that are smaller than the lower bound, chosen apriori, are made black. Correspondingly all 
pixel values larger than the chosen upper bound are set to 255, i.e., made white. The rest of the values that lie in between this 
minimum and maximum bound are spread linearly on a to 255 scale. 

The result of this attack on the watermarked image is shown in Fig. |15(a)| and that of the extracted watermark is shown in Fig. 
|15(b)| The enhancement value chosen in the present case is 0.75. We have also checked it for another value 1.25 and from the 
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FIG. 15. |(a)| Watermarked image after applying contrast enhancement and | (b) [ Extracted watermark. 




(a) (b) 



FIG. 16. (a) Dilated watermarked image and (b) Extracted watermark. 



PSNR and RMSE values given in the table, the quality of the embedded and the extracted watermark actually improves with a 
higher value of enhancement. 



K. Dilation 

Dilation in images can be thought of as a morphological operation. When it is applied to grayscale images, the bright regions 
that are surrounded by dark regions increase in size, whereas the dark regions having bright regions around it reduce in size. 
The effect is very prominent at places, in the image, where the intensity changes rapidly. Regions where the intensity is fairly 
uniform remain largely unaffected except at the edges. Thus the overall effect of this operation is to whiten the whole image. 
The results are shown in Figs. |16(a)| and |16(b)| 



L. Scaling 

Scaling operation changes the size of images. This can be of two varieties, namely, uniform and non-uniform. When the 
horizontal and vertical directions are scaled by the same factor it is called uniform scaling. In contrast to this, non-uniform 
scaling uses different scaling factors for the horizontal and vertical directions. This latter type of scaling changes the aspect ratio 
of the image unlike the former. The watermarked and extracted images after uniform scaling are shown in Figs. |17(a)| and |17(b)| 
respectively. The watermarked image. Fig. |17(c)[ and the extracted watermark. Fig. |17(d)| after non-uniform scaling show that 
our watermarking procedure performs quite well. 



M. Image Difference 

The absolute difference between the original image and the watermarked image, i.e. Fig. |l(a)| and Fig. |2(e)| is shown in Fig. 
[T8]reveals that the embedded watermark is completely invisible and we only see a corrupted texture of the original image. 



S.No. 


Type of Attack 


Software Used 


Dollar 


Pepper 


PSNR (dB) 


RMSE (E) 


PSNR (dB) 


RMSE (e) 


1. 


Tune Sharpen (Directional) 


Corel Graphics Suite 1 1 


23.56 


16.31 


13.19 


51.43 


2. 


Smoothing (50%) 


Corel Graphics Suite 1 1 


22.53 


18.36 


10.88 


67.14 


3 


Lowpass Blur (Radius 1) 


Corel Graphics Suite 11 


19.82 


25.09 


9.30 


80.47 


4 


Gaussian Noise 


MATLAB R2010a 


22.32 


18.81 


9.83 


75.71 


5. 


Gaussian Blur (Radius 1) 


Corel Graphics Suite 1 1 


19.18 


27.02 


8.95 


83.85 


6. 


Gamma correction 1.09 


MATLAB R2010a 


22.64 


18.13 


11.59 


61.81 


7. 


Gamma correction 0.95 


MATLAB R2010a 


24.92 


13.94 


15.17 


40.95 


8. 


Deinterlace (Even Une, Interpolation) 


Corel Graphics Suite 1 1 


18.87 


27.98 


11.58 


61.92 


9. 


Contrast Enhancement 0.75 


Corel Graphics Suite 1 1 


18.59 


28.92 


7.16 


103.04 


10. 


Contrast Enhancement 1.25 


Corel Graphics Suite 1 1 


24.00 


15.51 


8.30 


90.32 


11. 


Dilation 


GIMP 2.6.10 


13.86 


49.85 


8.28 


90.48 


12. 


lirosion 


GIMP 2.6. 10 


12.44 


58.73 


7.84 


95.26 


13. 


Local equalization (80, 80) 


Corel Graphics Suite 1 1 


11.90 


62.47 


6.38 


112.61 


14. 


Rotation 1 degree (Bilinear) 


MATLAB R2010a 


13.18 


53.89 


7.91 


94.48 


15. 


Additive Noise 


MATLAB R2010a 


18.56 


29.01 


13.37 


50.59 


16. 


Cropping 


MATLAB R2010a 


11.21 


67.66 


13.62 


48.98 


17. 


Median I'illcr (3x3) 


MATLAB R2()l()a 


18.86 


28.04 


9.8 


75.85 


18. 


Intensity Adjustment 


MATLAB R2010a 


20.86 


22.25 


9.52 


78.51 


19. 


Uniform Rescaling (512-256-512) 


Xn View 


19.17 


27.05 


9.06 


82.75 


20. 


Non-uniform Rescaling (512-(320x240)-512) 


Xn View 


5.78 


126.37 


9.10 


82.38 


21. 


JPEG 20% 


MATLAB R2010a 


20.70 


22.67 


9.67 


77.37 


22. 


JPEG 40% 


MATLAB R2010a 


21.98 


19.57 


10.47 


70.85 


23. 


JPEG 60% 


MATLAB R2010a 


22.76 


17.89 


11.08 


65.61 


24. 


JPEG 80% 


MATLAB R2010a 


23.63 


16.18 


12.66 


54.68 


25. 


JPEG2000 (Compression ratio 2) 


MATLAB R2010a 


24.26 


15.05 


20.23 


22.88 


26. 


JPEG2000 (Compression ratio 4) 


MATLAB R2010a 


23.92 


15.66 


13.55 


49.33 


27. 


JPEG2000 (Compression ratio 6) 


MATLAB R2010a 


23.26 


16.88 


11.50 


62.51 


28. 


JPEG2000 (Compression ratio 8) 


MATLAB R2010a 


22.64 


18.14 


10.50 


70.11 


29. 


JPEG2000 (Compression ratio 10) 


MATLAB R2010a 


22.07 


19.36 


10.01 


74.17 


30. 


JPEG2000 (Compression ratio 12) 


MATLAB R2010a 


21.51 


20.65 


9.68 


77.06 


31. 


JPEG2000 (Compression ratio 14) 


MATLAB R2010a 


21.09 


21.69 


9.42 


79.36 


32. 


JPEG2000 (Compression ratio 16) 


MATLAB R2010a 


20.72 


22.63 


9.19 


81.56 


33. 


JPEG2000 (Compression ratio 18) 


MATLAB R2010a 


20.43 


23.40 


9.05 


82.89 


34. 


JPEG2000 (Compression ratio 20) 


MATLAB R2010a 


20.19 


24.06 


8.93 


83.96 



TABLE 1. Table showing different types of attacks performed on the watermarked images. The PSNR and RMSE values of the watermarked as well as the extracted image are given. 
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(c) (d) 



FIG. 17. 1(a) I |(b)| Uniformly rescaled watermarked image and Extracted watermark [(c)] [(d)] Non-uniformly rescaled watermarked image and 
Extracted watermark. 




FIG. 18. Image difference between Fig. |l(a)| and Fig. |2(e)| 



VI. SUMMARY 



In this paper we have introduced a new watermarking scheme based on the Hilbert transform of digital images. Due to inherent 
ambiguities with the definitions of two-dimensional Hilbert transform, we apply the one-dimensional Hilbert transforms column- 
wise to any image matrix in which information is sought to be embedded. The Hilbert transformed digital image is defined 
through its matrix of phase and amplitude values. The main idea behind the proposed technique is that the phase of the analytic 
signal associated with typical digital images is generally a sparse matrix and hence offers large space for hiding information in 
a highly imperceptible manner. 

We have described embedding and extraction algorithms. The embedding is performed in the Hilbert transformed domain. 
The Hilbert transform of a watermark image is embedded in the Hilbert transformed domain of the host image, after suitable 
scaling by a factor X. An efficient algorithm for computing the optimal scaling factor X has been derived. The quality of the 
watermarked image and extracted watermark have been shown to be good when measured in terms of PSNR and RMSE. In 
addition, we have performed a large number of attacks to study the robustness of our algorithm. The results are summarised 
in Table The proposed method is shown to be robust against additive noise, cropping, gaussian noise, JPEG/JPEG 2000 
compression, median filtering, rotation, gamma correction, intensity adjustment, gaussian blur, contrast enhancement, dilation, 
rescaling and image difference. We would like to emphasise that calculating Hilbert transforms can be efficiently done on 
personal computers too. It is not computationally exacting. Thus, we believe that our algorithm can also be modified to suit 
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real-time applications in an efficient manner. 
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